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Abstract 

The sign problem is notorious in Monte Carlo simulations of lattice QCD with 
\Q . the finite density, lattice field theory (LFT) with a 9 term and quantum spin 

models. In this report, to deal with the sign problem, we apply the maximum 
1 entropy method (MEM) to LFT with the 9 term and investigate to what extent 

■ the MEM is applicable to this issue. Based on this study, we also make a brief 

comment about lattice QCD with the finite density in terms of the MEM. 

CO 

1 Introduction 

It is an important subject to reveal the phase structure of QCD in [i-T space, where T and \i are tempera- 
\ ture and quark chemical potential, respectively. This gives hints not only to understand the physics of the 

early universe and the neutron star, but also to analyze what happens in heavy ion collisions. The lattice 
simulation is one of the most reliable methods to comprehensively study the phase structure. However, 
Monte Carlo (MC) simulation based on the importance sampling method cannot directly apply to Lat- 
tice QCD at the finite density, because the fermion determinant with fi makes the Boltzmann weight 
complex. This is the notorious sign problem. Although various techniques to circumvent this problem 
have been proposed,! 1 1 the sign problem has not been solved yet. In this report, the maximum entropy 
method (MEM)[2 3] is introduced from a different viewpoint. By applying the MEM to lattice field 
theory (LFT) with a 9 term, where it also suffers from the sign problem, we investigate to what extent 
the MEM is applicable to this issue. Based on this study, we make a brief comment about lattice QCD at 



the finite density in terms of the MEM. 



X ' 2 Sign Problem in LFT with the 9 Term 

The partition function Z{9) in LFT with the 9 term can be calculated by Fourier-transforming the topo- 
logical charge distribution P(Q): 

where S(z, z) and Q(z, z) are the action and the topological charge as functions of lattice fields z and 
z, respectively. Note that P(Q) is calculated with a real positive Boltzmann weight. We call this the 
Fourier transform method (FTM). Although this method works well for small volumes, it breaks down 
for large volumes. This is because the error in P(Q) disturbs the behavior of the free energy density 
f(9) = — y log Z(9) (V is a volume). Figure [T]displays f(9) obtained from MC data of the CP 3 model 
with the fixed point action. The coupling (5 is fixed to 3.0 and various lattice sizes L are employed. 
The number of measurements is several millions for each case. Although f(9) for L < 38 behaves 
smoothly in the whole 9 region, f{9) for L = 50 and 56 cannot be properly calculated for 9 > 2.0. In 
the L = 56 case, f{9) becomes flat for 9 > 2.0. This is called flattening. In the L = 50 case, f{9) 
cannot be obtained for 9 > 2.0 due to negative values of Z{9). We also call it flattening, because the 
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Fig. 1: Free energy density f{6) obtained from the MC data of the CP 3 model. The coupling (3 is fixed to 3.0 and lattice sizes 
L are changed from 12 to 56. The number of measurements reaches several millions for each case. 



error in P(Q) causes this behavior in the same way as the L = 56 case. Flattening is originated from 
the sign problem. This is understood in the following way.|4]|5]] The MC data of P{Q) consists of the 
true value of P(Q), P{Q), and its error, AP(Q). When the error in P(Q) at Q = dominates, f{9) is 
approximated by f(9) ~ -y log[e _l//(9) + AP(0)]. Here, f(9) denotes the true value of f(6). Since 

f{9) is an increasing function of 9, e Vf ~^ ~ |AP(0)| could occur at 9 = 6 t and f(9) ~ - y log5P(0) 
for 9 > 9{. To overcome this problem requires the number of measurements proportional to e v . 



3 MEM 



The MEM is one of the parameter inference based on B ayes' theorem and derives a unique solution by 
utilizing data and our knowledge about the parameters. l2l l3l l6l In our MEM analysis. lf71 IH1 the inverse 
Fourier transform 

pw p -i8Q 

P(Q)= jj9—Z(9). (2) 

is used. The MEM involves to maximize the posterior probability prob(Z(9)\P(Q), I). Here, 
prob(^(c?)|i- > ((5), /) is the probability that Z{9) is realized when the MC data of {P(Q)} and informa- 
tion / are given. Information / represents our state of knowledge about Z{9) and Z{9) > is imposed. 
The probability is given by 



prob(£(0)|P(Q),I) ocexp 



\x 2 + aS 



= e 



W[Z] 



(3) 



where x 2 , a and S denote a standard x 2 -function, a real positive parameter and an entropy, respectively. 
Conventionally, the Shannon- Jaynes entropy 



S 



(19 



Z{9) - m{9) - Z{9) log 



Z{B) 
m{9) 



(4) 



is employed. A function m{9) is called default model and is chosen so as to be consistent with /. 
The most probable image Z{9) is obtained according to the following procedures. (1) To obtain the 
most probable image for a given a, Z^ a \9), by maximizing W[Z\. (2) To obtain the a-independent 
most probable image Z{9) by averaging Z^ a '(9) over a; Z(9) = J da P(a)Z^(9). The probability 
P(a) represents the posterior probability of a. (3) To estimate the error in Z{9) as the uncertainty 
of Z{9). The probability P(a) is given by P(a) oc g{a)e w ^ +A{ ~ a \ Here, W(a) = W[Z^], and 
A(a) represents contributions of fluctuations of Z(9) around Z^ a \9). The function g{a) is the prior 
probability of a. Conventionally, two types of g(a) are used: ^Lap(a) = const (Laplace's rule) and 
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ffjef (°0 = l/o (Jeffrey's rule). Information about a before obtaining data does not play the conclusive 
role in the derivation of Z(9). In the present study, the g (a) -dependence of Z(9) is estimated by the 
following quantity: 

|^Lap(#) - ^Jef 



A(0) = 



^Lap(^) 



(5) 



where Z\^ ap {9) and 2j e f (0) are the most probable images for Laplace's and Jeffrey's rules, respectively. 
4 Numerical Results 

We apply the MEM to the MC data with flattening as well as without flattening. The latter is the data 
for L = 38 (data A) and the former is those for L = 50 (data B). Here, two types of m{6) are used: (i) 
Gaussian function mc(S) = exp[— 7 lo ^. 10 6 )2 ], where a parameter 7 is changed over a wide range, and 
(ii) m{9) = Z{9) for smaller volumes. In case (ii), to analyze the data for L = L , Z{9) obtained by the 
MEM for smaller volumes are utilized as m{9). For Lq = 50, Z{9) for L = 24, 32 and 38 are used as 
m{9), which are denoted as m L / Lo (9) = m L / 50 (9). In this report, all results of the MEM with Laplace's 
rule are shown except for A(6>). In the analysis, the Newton method with quadruple precision is used. 

4.1 Non Flattening Case 

Figure 13 displays Z{9) for data A. The Gaussian defaults mo,{9) with 7 = 0.6 and 1.0 are used. The 
partition function Zf ou1 .(9) obtained by the FTM is also plotted. Both the results of the MEM have 
no m(#)-dependence and are in agreement with the result of the FTM. The error of Z{9), 5Z{9), are 
calculated according to the procedure (3). These errors are too small to be visible in Fig. |2] 
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Fig. 2: The most probable images in the non flattening case (L — 38). The result of the FTM is also plotted (x). 



4.2 Flattening Case 

Let us turn to data B. Unlike data A, much care is needed in the analysis.! 7 ] In order to properly evaluate 
Z(9) as the final image, we investigate (i) the statistical fluctuation of Z(9), (ii) g(a)-dependence of 
Z{9) and (iii) the relative error of Z{9). In (i), it is found that the statistical fluctuation of Z{9) becomes 
smaller with increasing the number of measurements and that Z(ff) with 20.0M/set is obtained with 
sufficiently small fluctuations except for near 9 = n. In (ii), we systematically investigate the g(a)- 
dependence of Z{9) by calculating A(#). The left panel of Fig. [3] displays A(#) at 9 = 2.60, as a 
representative. Here, the Gaussian defaults are used, where 7 is changed from 3.0 to 13.5. The value of 
A(#) is the smallest for 7 = 5.0 and becomes larger as the value of 7 deviates from 5.0. Similar results 
are obtained in the whole 9 region. This seems to indicate that rriQ,(9) with 7 = 5.0 is the most suitable 
as m{9) among the defaults which we have chosen. Keeping in mind that A(6>) includes an uncertainty 



4 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




Fourier 
y=3.0 
y=4.0 
Y=5.0 

m 24/50^) 
m 32/50 (° 
m 3S«n(o) 



mum* 



14 



0.5 



1 



1.5 



2.5 



3.5 



2 4 6 8 10 12 

Gaussian function with y 6 
Fig. 3: Values of A(0) at 9 — 2.60 (left panel). The Gaussian defaults are used. Values of \SZ(6)\/Z(9) for the selected 6 

default models (right panel). The value of \8Zf out (8)\/ Zf oui (8) is also plotted in the right panel. 



originated from 5Z{9), we impose a constraint that the final images should satisfy A(f5) < 0.2. Here, 
this value is chosen as a typical one of the uncertainty in A(0) coming from 5Z(9). Six images satisfy 
this constraint among those which we have obtained, and do not depend on m{9) up to 9 = 3.0. In (iii), 
we investigate how the MEM is applicable to our issue by calculating the relative error \5Z{9) \/Z(ff). 
Upon a constraint \5Z(9)\/ Z(9) < 0.3, the four most probable images Z{9) satisfy the constraint up 
to 9 = 3.0. This constraint is chosen from the fact that the error propagation of P(Q) starts to strongly 
affect the behavior of Zf oui (9) at \SZf oui (6)\/ Zf out (9) ~ 0.3 in the FTM (see the right panel of Fig. 
0). Here, this value realizes at smaller value of 9, 9 = 2.4. These results are displayed in the right panel 
of Fig. El 

In this analysis, we find that the four most probable images are obtained with reasonably small errors 
in a wide range of 9, which is displayed in the left panel of Fig. |4] As a comparison, Zp our (9) is also 
shown in the right panel. 
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Fig. 4: The most probable images for the selected m(9) (left panel). As a comparison, Zfout(9) obtained by the FTM is also 
displayed in the right panel. The total number of measurements is 30.0M/set for both the cases. 



5 Summary and Discussions 

In this report, to deal with the sign problem in LFT with the 9 term, we apply the MEM to the MC data 
of the CP 3 model. In non flattening case, all results of the MEM agree with the one of the FTM within 
the errors. In the flattening case, obtained images depend on m{9). By investigating whether they are 
adequate images, we have found that the MEM allows us to calculate Z{9) with small errors for large 9 
region. For the details, see Ref. [9| 

Finally, let us make a brief comment about lattice QCD with the finite density in terms of the MEM. 
In lattice QCD with the finite chemical potential, MC simulation cannot be directly performed due to 
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the complex phase of the fermion determinant. There are various techniques to avoid the sign problem 
and we concentrate on the canonical ensemble approach. 1 101 ITTI [T2l [T3l By the fugacity expansion, 
Z(V, T, fj,) is written as 

Z(V,T^) = J2z(V,T,n)(e^ T ) n , (6) 

n 

where n is the total quark number. Taking /x = icj)T, where ^ is a real, Z(V, T, \i = icj)T) is free 
from the sign problem and Z(V,T,n = i(j)T), in principle, can be calculated with MC simulation. In 
this case, Z(V,T,fi = ifiT) = J2n %(V, T, n)e in ^. Comparing it with Eq. ©, we see the following 
correspondence: 

{P(Q) «-► Z(V,T,fi = i(j)T), e- ieQ /2ir «-> e^ n , Z(9) <-» Z(V,T,n)}. (7) 
It may be worthwhile to study the theory in terms of the MEM. 
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